Sub-natural linewidths in excited-state spectroscopy 
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We investigate, theoretically and experimentally, absorption on an excited-state atomic transi- 
tion in a thermal vapor where the lower state is coherently pumped. We find that the transition 
linewidth can be sub-natural, i.e. less than the combined linewidth of the lower and upper state. 
For the specific case of the 6P3/2 — ^ 7Si/2 transition in room temperature cesium vapor, we measure 
a minimum linewidth of 6.6 MHz compared with the natural width of 8.5 MHz. Using perturba- 
tion techniques, an expression for the complex susceptibility is obtained which provides excellent 
agreement with the measured spectra. 
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I. INTRODUCTION 
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Spectroscopy of excited state transitions is of growing 
interest for a variety of applications including the search 
for stable frequency references [TJ [5] , state lifetime mea- 
surement [3] , optical filtering ^ , frequency up-conversion 
[5], multi-photon laser coohng [6], as well as Rydberg 
gases [3 IH] and their application to electro-optics [5HTT] 
and non-linear optics fT51. Excited state spectroscopy 
can be achieved without significant transfer of population 
out of the ground state using electromagnetically induced 
transparency (EIT) [T3| in the ladder configuration [H] . 
In conventional EIT, the excited state transition is driven 
by a strong coupling laser creating a transparency win- 
dow which is then detected using a weak probe on the 
ground state transition. In thermal vapors, ladder EIT 
is only possible when the lower transition is probed [15] 
and the probe wavelength is greater than the coupling 
wavelength [TB] . Alternatively, on strong transitions such 
as the infra-red transitions from excited states in alkali 
atoms, one can probe directly on the excited state tran- 
sition and detect absorption or fiuorescence [17]. 

In this work we develop the theory for the complex 
susceptibility of excited state transitions and compare 
the results to experimental observations. We consider 
a transition between two states, |1) and |2) shown in 
Fig. 1 (a). Neglecting the effects of Doppler broaden- 
ing, it is expected that the lineshape of the absorption 
is a Lorentzian whose FWHM is given by the sum of 
the natural linewidths of states |1) and |2), Fi + F2, as 
shown in Fig. 1(b). We show that if the population of 
state |1) is coherently pumped from another state with 
a weak coupling field, Fig. 1(c), the absorption lineshape 
remains a Lorentzian but its FWHM is solely determined 
by the natural linewidth of the upper state, F2 as shown 
in Fig. 1(d) because the lower state is effectively stable. 

The paper is organized as follows. In Section 2 we 
develop the equations of motion for the system and de- 
rive an expression for the susceptibility. In Section 3 
we extend the analysis to include the effects of Doppler 
broadening. In Section 4 we compare our derived expres- 
sions to experimental data in the limit of weak pumping. 
In Section 5 we compare theory and experiment in the 





FIG. 1: (Color online) (a) Schematic of the energy levels of 
an atom whose lower state 1 1) is incoherently populated. The 
natural decay of states |1) and |2) are Fi and F2, respectively. 
The absorption lineshape is probed by the probe field of Rabi 
frequency Q.p. (b) The FWHM of the lineshape is Fi -I- 
(c) Schematic of the energy levels of a transition whose lower 
state 1 1) is coherently populated by the coupling field, of Rabi 
frequency Q,c, from a state |0). In this case the FWHM of the 
absorption lineshape is r2 when is small, as shown in (d). 



regime of strong coupling before concluding in Section 6. 
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II. EXCITED STATE TRANSITION WITH 
COHERENT PUMPING OF LOWER STATE 



tions of motion for the density matrix elements, Py, are 



Equations of motion and the steady state 
solutions 



Consider a transition between 2 states as shown in 
Fig. [ijc), a lower state, and an upper state, |2), 
with the associated eigenenergies of huji and h,uj2, re- 
spectively. Initially, states |1) and |2) are not populated. 
To populate state |1), the system is coherently pumped 
by the resonant coupling field with Rabi frequency of 
ric from the stable eigenstate |0) whose eigenenergy is 
fiujo and wq < wi < uj2- The transmission (or absorp- 
tion) lineshape of the transition is probed by scanning 
the frequency of the probe field whose Rabi frequency is 
rip. Applying the rotating- wave approximation and the 
slowly- varying variables transformation, the Hamiltonian 
of the system is given by IK = 3-Cq + J£i, where 
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with Ap — LUp — (uj2 ~ uji) , Ac 



(cji — loq) , and 



Ar = Ap -I- Ac . Here Ap^(c) is the detuning of probe 
(coupling) laser, Wp_(c) is the angular frequency of probe 
(coupling) laser, Ar is the two-photon Raman detuning, 
and h.c. is the hermitian conjugate. The first term of the 
total Hamiltonian, JCg, represents the field- free atomic 
system, whereas the second term of the Hamiltonian, !Ki, 
describes the interaction with both probe and coupling 
fields. 

Using standard semiclassical methods [T^, the equa- 
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- ^Pi2, (2f) 
where we define effective linewidths 7' = ri/2-l-7c, 7" — 
(Ti + r2)/2 + 7p , and 7"' = r2/2 + 7p + 7c and Ti 
and r2 are the natural linewidths of the states |1) and 
1 2) respectively. In addition to spontaneous decay, we 
include a dephasing due to the linewidth of the probe 
and coupling fields of 7p and 7c, respectively. Solving 
Eqs. ([2| (with pij — 0) together with the constraint 
Poo + Pll + P22 = 1 using a perturbation technique (see 
Appendix [A| , the steady state solutions of the density 
matrix pij are given by. 
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In the weak excitation limit of flp <ti 7", P22 — 0. We can 
therefore assume that no population is lost from the sys- 
tem via other decay channels from state |2), although ad- 
ditional decay channels may contribute to the linewidth, 

r2. 

We consider pi2 as it determines the complex suscep- 
tibility of the system. It is clear from Eq. (3d) that the 



number of atoms pumped into state |1) strongly affects 
the magnitude of pi2, manifest as the multiplication fac- 
tor piirip/2. In our case, the population of state |1) is 
resonantly pumped from state |0), i.e., Ac ~ 0, and 7c 
is much less than 7". Thus, the second term in the first 
square bracket approaches unity within this approxima- 
tion. 
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B. The complex susceptibility of the system 

The complex susceptibility of the system at the probe 
frequency is obtained by comparing the polarization ob- 
tained from classical electrodynamics with that calcu- 
lated using a density matrix treatment [13l [14] . The ex- 
pression of the complex susceptibility is then given by 

m 



X = 



-P2i 



(4) 



where Ji is the atomic density and d2i is the dipole ma- 
trix element for probe transition. For |x| < 1 the real 
part and imaginary part of the susceptibility are respec- 
tively proportional to the refractive index tir and the 
absorption coefficient ex. via the relations 
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where is the wavevector of the probe field. Thus the 
complex susceptibility of the system is 
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Here J^i = ^Pii, the atomic density of state Accord- 
ing to Eq. (6) , the complex susceptibility is proportional 
to the number of the atoms in state As state |1) 
is coherently pumped from state |0), the multiplication 
factor has a Lorentzian profile as a function of A^.. This 
implies that to get a large susceptibility, one needs to 
resonantly pump population to state The term in 
the square bracket is similar to the result obtained for 
EIT in the three-level cascade system [T3]. It indicates 
that, for a finite value of flc, the system will become 
transparent when the probe field is scanned across the 
resonance at Ap = due to Autler-Townes splitting [20] . 
When state |1) is weakly pumped by the coupling field, 
i.e., ill ^ r2(r2 + Ti), the term in the square bracket 
of Eq. ^ can be expanded using a Taylor expansion. 
Neglecting the higher order terms in flc, the complex 
susceptibility reduces to 
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This complex susceptibility is similar to that of a 
two-level system, except for the multiplication factor. 
The transmission lineshape of the system is simply a 
Lorentzian centered at —Ac with FWHM 7"'. For a suf- 
ficiently small laser linewidth compared to the ap- 
proximation of the FWHM is solely determined by the 
linewidth of the excited state, irrespective of the 
linewidth of state Fi, i.e.. 
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This result is different to the case in which state |1) is 
incoherently populated. In such a case, the FWHM of 
the transmission lineshape is determined by the sum of 
the linewidths from both the lower state and the upper 
state, i.e., Fi -I- F2 [21,. 

For experiments in room temperature vapors, the 
Doppler effect must be included into the model and this 
topic will be discussed in the next section. 



III. EFFECT OF DOPPLER BROADENING 

Each atomic velocity class in a thermal vapor experi- 
ences a different laser detuning Ap and Ac due to the 
Doppler effect. To obtain the velocity-dependent com- 
plex susceptibility, we make changes to Ap, Ac and 3Nf 
with the following substitutions [14j : 
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where fcp(c) is the wavevector of the probe (coupling) 

field, u — ^Ik-QT jm is the most probable speed of the 
atoms at a given temperature T and m is the mass of an 
atom. Substituting Eqs. ([9]) into Eq. the value of 
the complex susceptibility of a particular velocity class v 
is then given by. 



X(i;)dw 



heo^k^ik, - kp)u^ 2Fi (z + + 



z - zo + 



17cV4 



(/cc - kp)kpu'^{z - zi) 
with the change of variable z = v/u, and 



7 



(^C ^p)'^ 




y2 _^ 



Zl 



ike kp)u 
A^ 

fccU ' 

- 17 , 



A„ + i7" 



, (10) 

(11a) 

(lib) 
(11c) 

(lid) 
(lie) 
(llf) 



The total susceptibility is obtained by integrating Eq. 
( 10 1 over all velocity classes. The full result of the in- 



tegration is discussed in Appendix B. We consider the 
case in which the coupling Rabi frequency, f^c is suffi- 
ciently weak that the EIT-like third term in the square 
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bracket of Eq. ( 10 ) is neglected. In this case the complex 
usceptibi 
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susceptibility, xd, becomes 
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From Eq. ( 12 1, the total complex susceptibility is simply 



the convolution between a Lorentzian of width 7 (term in 
square bracket describing the transition from lower state 
|1) to upper state |2)) and a product of a Lorentzian of 
width a and a Gaussian (term in curly bracket describing 
the atomic velocity distribution of lower state |1)). 

The result of the integration in Eq. ( 12 ) involves the 
Faddeva function (the exact result of the integra- 
tion is described in Appendix C). However, the integra- 
tion can be simplified by replacing the product between 
a Gaussian and a Lorentzian with a Lorentzian, i.e., 



(13) 



This approximation is valid since, at room temperature, 
the width of the Lorentzian a is much smaller than the 
width of the Gaussian. Hence the product of the Gaus- 
sian and the Lorentzian approximately vanishes when 
\z\ > a and we can approximate the product by the 
Lorentzian of width a. In other words, the range of 
velocity classes involved in the integration around the 
position where the Lorentzian of width a is centered 
is much smaller than the most probable speed of the 
atoms, V <^ u. Hence, we can expand the Gaussian 
around the position where the Lorentzian is centered, 
i.e., exp(— z^) « exp(— /J^). 

Using this approximation, the total susceptibility is 
simply given by 
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It is clear from Eqs. (14) and (15) that the absorption 
lineshape remains Lorentzian with the FWHM of Ffwhm 
even when the Doppler effect is included. It is worth not- 
ing that the total susceptibility in this case is different 
from the total susceptibility calculated for the case of 
incoherent pumping. The total susceptibility for inco- 
herent pumping is the convolution between a Lorentzian 
and a Gaussian, resulting in a Voigt profile [23]. 

In the limit where ^c/^i ^ 1, the linewidth of the 
absorption profile is simply. 
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(neglecting 7p and 7c). It contains the sum of two terms: 
the first term is the linewidth of the absorption lineshape 
in the case in which the Doppler effect is neglected and 
the latter is the linewidth of the lower state scaled by 
the ratio obtained from the wavevectors. Physically, the 
second term originates from the fact that the atoms are 
velocity-selected by the Doppler effect for the atom-field 
interaction. Only atoms whose velocities are between 
—Ti/2kc and ri/2fcc are coherently pumped into the 
lower state when the coupling field is on resonance. Since 
Ti/kc is very small compared to the width of the Doppler 
broadening, all atoms pumped into state 1 1) have approx- 
imately the same velocity distribution, i.e., the distribu- 
tion is independent of velocity, and given by 



fiv) 
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However, the distribution of the atoms is also determined 
by the Lorentzian of width Fi/Zcc- Thus the final veloc- 
ity distribution of the pumped atoms is a Lorentzian of 
width Fi/Zcc and the height is scaled by Eq. (17). In the 
two-photon interaction process, the width of Ti/kc in ve- 
locity space is equivalent to the width of (fcc — A:p)ri/fcc 
in frequency space. Hence the total linewidth of the final 
absorption lineshape is the sum of (fcc — fcp)ri/A:c with 
the unaffected linewidth, 



IV. COMPARISON BETWEEN THEORY AND 
EXPERIMENTAL RESULTS 

To test the result derived above we used the experi- 
mental setup described by Carr et al. [21] • The experi- 
ment was performed in a 7.5 cm vapor cell containing Cs 
at room temperature. The 1470 nm weak probe beam 
(with horizontal linear polarization) and the 852 nm co- 
axial, counter-propagating coupling beam (with horizon- 
tal linear polarization) are applied along the vapor cell 
axis. The probe and coupling beams have l/e^ radii of 
1.2 mm and 1.6 mm, respectively. The coupling beam 
was stabilised to the 6Si/2,i^ = 4 — > 6P3/2,F' — 5 
transition while the probe beam was scanned across the 
6P3/2,-F' = 5 7Si/2,F" = 4 transition. The scan is 
calibrated using a wavemeter to better than 1% accuracy. 
The transmission signal measured from the experiment is 
shown as the solid black line in Fig. [2ja). 

To model the transmission lineshape, the complex sus- 
ceptibility is calculated for each magnetic sublevel and 
the total complex susceptibility is the average of all com- 
plex susceptibilities over all possible magnetic sublevels. 
This is given by 



XTOT(Ap) = ^ ^ xr(Ap), 



(18) 



m j='——4 



where ^ (Ap) is the complex susceptibility correspond- 
ing to the niF magnetic sublevel of 681/2,-^ = 4 state. 
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FIG. 2: (Color online) (a) Comparison between the observed 
spectra, shown by a black solid line, for a relatively small 
Ql — 0.6 MHz, and the theoretical model. The grey (blue) 
solid line is the theoretical prediction, taking into account the 
absorption of the coupling field across the vapor cell, (b) The 
residual plot between the observed data and the theoretical 
model. 
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FIG. 3: (Color onhne) The values of the FWHM are plot- 
ted as a function of 0,1. The dashed curve is the theoretical 
prediction calculated using the weak coupling model, whereas 
the solid (red) curve is the theoretical prediction calculated 
using the full model. As expected the data are in agreement 
with the weak coupling model when is relatively small. 
However, for a large Ql, the weak coupling model fails to pre- 
dict the FWHM. This difference at large Ql can be recovered 
when using the full model. The combined lower and upper 
state linewidth Fi -I- F2 is shown by the dotted line. 



Xo" is calculated using Eq. (14), where the coupling 
Rabi frequency is sufficiently weak. The factor of 1/16 
in the equation accounts for the fact that the initial pop- 
ulation is evenly distributed among the magnetic sub- 
levels of 6Si/2,-F — (3,4). The coupling Rabi frequency 
of the transition 681/2,^ = 4 6P3/2,F' ~ 5 and 
the dipole matrix element of the transition 6P3/2 7^" = 
5 — > 7Si/2, F" = 4 corresponding to each magnetic sub- 
level, fi™^ and , are given by 

= 5.63eao x ( J^, J ) = (19b) 

where the reduced dipole matrix element of the transi- 
tion 6Si/2,-F = 4 — 7'6P3/2,-F" = 5 is absorbed into fiJJ, 
i.e., fij; = ei?c(6P3/2||?'||6Si/2)/?i, ao is Bohr radius, and, 
mp and mp' are the magnetic sublevels of the 6S1/2 and 
6P3/2, respectively. 

The comparison between the experimental data and 
the theoretical model is shown in Fig. 2(a). The theoret- 
ical transmission shown as the solid grey (blue) line was 
calculated using Eqs. (5b) and (18) for a known tem- 
perature of T = 22 °C. The fit parameters fll/2'K = 0.6 
MHz, 7p/27r = 0.2 MHz are in good agreement with the 
experimental values taking into account the variation of 
the coupling intensity across the vapor cell. The resid- 
ual plot in Fig. 2(b) shows that the theoretical model is 



in good agreement with the experimental data. Fig. 3 
shows the transmission lineshape FWHM as a function 
ill- The dashed line shows the FWHM as a function of 
riJJ for the weak coupling field approximation. The de- 
pendence on fll shows good agreement when fiJJ is small. 
However when f2JJ is large, the disagreement between the 
model and the experimental data increases. The agree- 
ment between the theory and the experiment can be re- 
covered when the complete solution (see Appendix B) is 
used to calculate the FWHM as a function of fiJJ (shown 
as the solid line). Note that both curves approach the 
same value of 27r x 6.6 MHz when ft^ approaches zero. 
This is sub-natural and less than the combined linewidth 
Fi + F25.2 + 3.3 = 8.5 MHz |15] (shown as the dot dash 
horizontal line). The difference between the model and 
the experiment at large f2JJ arises from Antler- Townes 
splitting [5D] and is described by the EIT-like term in 
Eq. (6). The comparison between experiment and the- 
ory for large il^ will be the topic of the next section. 



V. COMPARISON BETWEEN EXPERIMENT 
AND THEORY FOR LARGE Ql 

We now apply the same method to model the transmis- 
sion lineshape, except that we are no longer in the weak 
pumping regime and the third term in Eq. (10) can no 
longer be neglected. Eq. (B8) is now used to calculate 
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FIG. 4: (Color online) (a) Comparison between the observed 
spectra, shown by a black solid line, for Q.]./2tt = 3.0 MHz, 
and the theoretical model. The grey (blue) dotted line is 
the theoretical prediction calculated using the weak coupling 
model, whereas the grey (red) solid line is the theoretical pre- 
diction calculated using the full model, (b) The residual plot 
between the observed data and the theoretical model calcu- 
lated using the full model. 
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FIG. 5: (Color online) (a) Comparison between the observed 
spectra, shown by a black solid line, for f2c/27r = 15 MHz, 
and the theoretical model. The grey (red) solid curve is the 
theoretical prediction calculated using the full model, (b) The 
residual plot between the observed data and the theoretical 
model. 



VI. CONCLUSION 



the complex susceptibility. 

Fig. 4(a) shows the comparison between the observed 
transmission lineshape, shown by a black solid line, and 
the theoretical transmission lineshapes calculated using 
both the weak coupling approximation, shown by a grey 
(blue) dashed line, and the full model, shown by a grey 
(red) solid line. The theoretical curves were calculated 
with = 3.0 MHz. The transmission calculated 

using the full model shows good agreement with observed 
data, resulting in the small residual as shown in Fig. [4jb). 
It is observed that the weak coupling approximation is 
in good agreement with the observed data, except at the 
region around the resonance. 

Fig. 5(a) shows the transmission lineshape, shown by 
a black solid line, in a case where ilJJ is large enough to 
Autler-Townes split the absorption resonance. The value 
of r^c/^'"' in this case is 15 MHz. The grey (red) solid 
line is the theoretical prediction calculated using the full 
model in the region |Ap/27r| < 20 MHz. Both the the- 
oretical prediction and the observed data are in good 
agreement around resonance. Fig. 5(b) shows the resid- 
ual between the theoretical prediction and the observed 
data. 



We have developed the theory of absorption lineshapes 
on excited state transitions where the lower state is coher- 
ently populated. We show that for an atom at rest, in the 
limit of weak pumping, the lineshape is a Lorentzian and 
the linewidth of the transition reduces to the linewidth 
of the upper state. Including the effect of Doppler broad- 
ening the linewidth is still subnatural and we find that 
the predicted lineshape is in very good agreement with 
experimental data over a wide range of coupling field pa- 
rameters. 
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Appendix A: Steady state solutions by perturbation 
technique 

Since the probe field is sufficiently weak (it will be 
shown later that the weak probe condition is fulfilled 
when r2p/7' <C 1), one can consider the expansion of 
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the density matrix, pij, in the power of Jlp, namely, 

P^, = + P^'^^V + Pf^^l + pS^ + • ■ • , (Al) 



Similarly, the set of equations corresponding to the n^^ 
power of fip (for n> 1) is given by 



where p^"^ is the n}^ order correction of the expansion of 



Pij. To solve for pij, we substitute Eq. (Al) into Eqs. 
(|2|), equate the terms of the same power inilp, and then 

solve for p-J-* from n = to all n. [26] . 

Applying this technique to Eqs. ([2]), the set of equa- 
tions corresponding to the zeroth power of ilp is given by 
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We find that the zeroth order corrections of pij vanish, 
except Poo\ p^^\ ^'^^ P^m- The expression for pp^-* and 
p\{ are given by 
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+ Al7e(pi^)-p[,i;))=0, (A4c) 



Pr'^-pir'M=0, (A4d) 
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Using Eqs. (|A3|) and Eqs. (|A4|), all of p^^^ again vanish. 



except P02'' and ; whose expressions are given by, 
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("2) ("2") 

It can be shown that p]^2 = P02 = 0. The solutions of 
the coherence po2 and pi2, to second order, are then 



P02 — f^pPQ2 J 
P12 = f^pPl2' • 



(A6a) 
(A6b) 



Thus far, we have assumed only that the probe Rabi 
frequency is sufficiently weak without quantifying this 
condition. To quantitatively determine the weak probe 
condition, let us consider the steady state of poi as this 
quantity is strongly related to pn and hence also to pi2. 
Using Eq. ( 2d I with the substitution poo = 1 — pn — P22, 



the expression of the steady state of poi is given by. 
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The first approximation is to neglect the contribution 
to poi from P22 as P22 vanishes up to the second order 
correction. This approximation is justified as the probe 
laser is weak and the upper state population is negligi- 
ble. Thereafter we consider the last term which contains 
the product between (r2p/2)/(iAc -I- 7'/2) and po2. The 
first term is of order ^^j^' at resonance Ac = 0. The 
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secon d term po2 is also of order 0.^/^' at resonance (see 
Eq. (A5a) and Eq. (7a)). Thus the product of these 
two terms is of order (rip/Fi)^, which can be neglected 
if rip/7' is much less than one. It follows that these ap- 
proximations lead to the same results for pij given by 
Eqs. (8). Hence, the weak probe condition is valid when 

f^p/y « 1- 



Appendix B: Complete solution of the Complex 
susceptibility of the system 

The complete solution of the complex susceptibility 
can be found by evaluating the integral 



(z + /3)2+ct2 

z - zo + 



(fee - kp)kpu'^{z - zi) 



(Bl) 



To evaluate the integral, one re-writes the integrand, us- 
ing partial fractions, as 



(z + /3)2+a2 



»g/4 

(fcc - kp)kpu'^iz - zi) 



- (t)+) 



[(/3 + 0_)2 + a2](0+-0_)z-(/._ 

i(zi + /3 + icr) c-^' 
2ct(/J + 0+ + i(T)(/3 + 0_ + icr) z + /3 + i(T 



i(zi + /3 — icr) 



2(j{l3 4- 0+ - icr) (/3 + 0_ - icr) z + /3 - icr 



(B2) 



where 



0± =2*-^° ^ 



(B3) 



The integrand has four poles in the complex plane, i.e., 
at (^+, 0_, — /3 — icr and — /3 -I- icr. Clearly, this complex 
integral reduces to the integral of the form 



-dz 



(B4) 



where Zp is the pole in the complex plane. The solution 
of the integration is given by 



-dz = isnW^sZp) 



(B5) 



where s = Sgn[Im(zp)] and Sgn is known as the sign 
function and its value is -t-1 when the argument is positive 



and —1 when the argument is negative. W{z) is known 
as Faddeeva function and it is defined as 



W{z) — e ^ erfc(— iz) . 
Hence the complex susceptibility is given by 



(B6) 



XD = 



h£Qk^{kc — kp)u'^ 2Ti 

(zi - (f>+) 



-4>-) 



{zi-4>-) 



i(zi + /3 + icr) 



2cr(/3 + (f)++ ia){l3 + 0_ + icr) 
i(zi +/3-ia) 



-s_VF(s_0_) 



W{~l3~ia) 



2a{l3 + (t)+~ia){(5 + (t)_-ia) 



— W{-l3 + ia) 



(B7) 



where s+ = Sgn[lm(0_|_)] and s_ — Sgn[Im((/)_)] 

One can use the same approximation as discussed in 
the article to approximate Eq. (B7) and the approxi- 
mated expression of Eq. (B7) is given by 



XD 



-e X 



TiEQ^klik^ - kp)u^ 2Ti 
tt{/3 + zi)(0_ - (/)+) + is+TTa{(f>+ + 2zi) 
a{(j)+ - + 0+ + is+(j){(3 + </>_- is+a) 



(B8) 



Appendix C: The exact result of Eq. (13) 

To evaluate the integral in Eq. (13), the integrand can 
be re-written using partial fractions as 



72 



(z^ +a^){z + ^ + i'y) cr(cr - 7 -f iO(^ + icr) 

e-"V2 

cr(CT + 7 - i^)(z - icr) 



[a2 + (e + i7)2](z + e + i7) 



(CI) 



Using Eq. (B5) and re-arranging the expression, the com- 
plex susceptibility is given by 



XD(Ap) 



1 



4heok^{kc - fcpu3) + + i-fY 



2 ^ZC] 2 

e"^oerfc(-izo) H e'^ erfc(cr) 

cr 



, (C2) 



where zq = ^ -F 17. 
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